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We relate a large class of classical spin models, including the inhomogeneous Ising, Potts, and clock models 
of q-state spins on arbitrary graphs, to problems in quantum physics. More precisely, we show how to ex- 
press partition functions as inner products between certain quantum stabilizer states and product states. This 
connection allows us to use powerful techniques developed in quantum information theory, such as the stabi- 
lizer formalism and classical simulation techniques, to gain general insights into these models in a unified way. 
We recover and generalize several symmetries and high-low temperature dualities, and we provide an efficient 
classical evaluation of partition functions for all interaction graphs with a bounded tree-width. 
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1. — Introduction. Classical spin models, such as the 
Ising, Potts and clock models [ 1 ], are widely studied in statis- 
tical physics. Despite of their simplicity, these models show a 
number of highly non-trivial features, which is e.g. apparent 
in their rich phase structure and critical behavior. These mod- 
els are considered in the context of (anti)ferromagnetism and 
are related to a number of central problems in lattice statis- 
tics. A quantity of fundamental interest is the partition func- 
tion, from which other physically relevant quantities can be 
derived. General symmetry and duality relations, as well as 
the efficient evaluation of the partition function are important 
problems in this context. In particular, it is known that for gen- 
eral interaction patterns the evaluation of the partition function 
of e.g. the Potts model is NP-hard |2], rendering the problem 
highly non-trivial. 

In this letter, we relate general classical spin models of q- 
state spins that interact pairwise according to some arbitrary 
interaction pattern, to problems in quantum theory, thereby 
obtaining a novel approach to tackle these problems. More 
precisely, we establish a correspondence between, one the one 
hand, classical partition functions, and, on the other hand, 
overlaps between certain quantum states — known as stabi- 
lizer states 01 01 — and complete product states of N q- 
dimensional quantum systems; here N is the number of pair- 
wise interactions. The stabilizer states encode the interaction 
pattern of the model, while the product states encode the de- 
tails of the interaction. This correspondence allows us to use 
a whole body of powerful results obtained in the context of 
quantum information theory regarding the classical descrip- 
tion of such quantum states, including the stabilizer formal- 
ism_0,0,l3l and the classical simulation of quantum systems 
(a, [J. Conversely, the vast knowledge on classical spin mod- 
els yields insights in the possible use of the corresponding 
states in quantum information tasks. 

We emphasize that the current approach allows one to in- 
vestigate classical spin models by virtually going to quan- 
tum systems — this is not be to confused with the solution of 
classical problems by implementing quantum algorithms on 
a quantum computer. Despite of the fact that the formula- 
tion of the problem in terms of quantum theory appears to 



increase the complexity and the number of involved parame- 
ters, this approach, as we will show, leads to a fruitful way 
to gain insights in spin models. For example, recently estab- 
lished classical algorithms o [J that allow one to describe 
and locally manipulate certain quantum systems in an efficient 
way, can be utilized to efficiently evaluate the partition func- 
tions of these models for all interaction patterns correspond- 
ing to graphs with a (logarithmically) bounded tree-width 1 8]. 
Furthermore, our approach leads to a simple way of recover- 
ing and extending general symmetry and duality relations of 
the above models in a unified way. For instance, using the 
fact that the corresponding states are stabilizer states, one can 
relate the stabilizer groups of the states to symmetry groups 
of the corresponding partition functions. In the following we 
will first describe the general class of classical spin models 
we consider, and then formulate the problem in a quantum 
language and apply techniques established in the context of 
quantum mechanics to gain insight in the classical models. 

2. — Classical spin models. We will consider a general 
class of classical spin models on arbitrary graphs. We con- 
sider n classical spins {s;} that can assume q possible states, 
Si S {0, . . . , q — 1}, and which interact pairwise according 
to an interaction pattern specified by a graph. More precisely, 
let G = (V, E) be a (connected) graph with vertex set V as- 
sociated to the spins, and edge set E. Further, let G" be an 
oriented version of G, where, for every edge e, one end-vertex 
v + is assigned to be the head of e , and the other end-vertex 
v~ is the tail of e. The Hamiltonian of the system is given by 
the general expression 

Here every h e is a local Hamiltonian defined on the edge 
e with the only restriction that the interaction strength be- 
tween the spins on the endpoints w+ and v~ is a function of 
the difference between s + and s - modulo q, denoted by 
l s n+ — s u~ \q- Thus, a general inhomogeneous model on an 
arbitrary graph is obtained. 

Note that this general model specializes to known and 
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widely studied statistical mechanics models such as the stan- 
dard Potts model, the Ising model and the clock model 1 1]. 
For example, the q-state Potts model (which specializes to 
the Ising model when q — 2) is obtained by taking h e (j) := 
— J e S(j, 0), for every j E {0, ...,</— 1} and for every e 6 E, 
where the J e are real coefficients. 

The central quantity in this context is the partition function, 

Z G (q,a,{h e }):=J2e- m{si}) , (2) 

where (3 — (ksT)^ 1 , with ks the Boltzmann constant and T 
the temperature, from which other relevant physical quantities 
such as e.g. free energy can be derived. 

3. — Quantum formulation. We now reformulate the par- 
tition function Zq in terms of a quantum physics language. 
To this aim, we will associate with the oriented graph G a a 
quantum state of N q-dimensional systems, where N := \E\ 
is the total number of edges of G, in the following way. Con- 
sider the incidence matrix B a of the oriented graph G a 
This is the matrix with entries and ±1 where the rows are 
indexed by the vertices of G, and the columns by the edges, 
such that (-B CT ) a .e = 1 if the vertex a is the head of the edge 
e, {B <7 ) a ,e = -1 if a is the tail of e, and (B cr ) a . e = oth- 
erwise. Let C G (q,a) be the row space of B a when arith- 
metic is performed over 1 q . This space is a linear subspace of 
:= Q) eeE 1 q - We now associate to C G (q, cr) the quantum 
state 

Kfe>:= £ \v), (3) 

where the \v) are basis states of (C 9 )®^ obtained by taking 
tensor products of a local basis {|0), . . . , \q — 1)}. Note that 
we have omitted the dependence of \ip G ) on q and cr to keep 
notation simple. 

An interesting property is that the state \4>g) is a stabilizer 
state, i.e., it is the unique joint eigenstate of q N commuting 
Pauli operators flj, ll ill . More specifically, defining the opera- 
tors X and Z by 

X\j) = \j + lmodq), and Z\j) = e^\j), (4) 

for every j = 0, . . . , q — 1, one finds that \ipa) is invariant 
under the q N local operations 

X(u)Z(v) :=(g)X u °Z v °, (5) 

eeE 

for every u 6 C G (q, cr) and v G C G (q, cr) . It is important to 
remark that the above construction of associating a quantum 
stabilizer state \ip G ) to a graph G is entirely different from 
the standard way of defining a graph state, usually denoted 
by \G) H. While for \tp G ) the ed ges of G are identified with 
quantum systems and the incidence matrix of G is used, yield- 
ing an TV -particle state (where N = in the construction 



of \G) the vertices are identified with quantum systems and 
the adjacency matrix is used, leading to a state on n = \V\ 
particles. 

We can now relate the partition function Zq to the state 
\i/jg), and obtain the first main result of this letter: The par- 
tition function Zc(q, cr, {h e }) on an arbitrary graph G with 
N edges can be written as the overlap between a quantum 
stabilizer state on N q- dimensional systems, and a complete 
product state. More precisely, we find 

Z G (q,a,{h e })=q-(^ G \ ((g) |a e ) J , (6) 

\e£E ) 

where \a e ) :— J^jZo e ^ l3h '^' 1 \j) is a c/-dimensional vector 
associated to the edge e, specifying the interaction 1 12]. 

Note that, whereas the states \ip G ) encode interaction pat- 
terns, the product states \a e ) encode the details of the interac- 
tions. For example, in the case of the Potts model one has 

|a e >=e-^|0) + (|l) + -.. + | ? -l)). (7) 

Examples of states \tp G ) are obtained as follows. First, if 
G is a tree graph — for example, a ID chain with open bound- 
ary conditions — it is easy to show that C G (q, cr) is the com- 
plete space (independent of the orientation cr). One then 
finds that the corresponding quantum state is a product state, 

namely \ip G ) = \+)® N , where |+) := |0) H h \q - 1). 

Hence, the partition function Z G on trees can easily be eval- 
uated, and one has Z G ~ JT (+|a e ). Second, a ID chain 
with periodic boundary conditions corresponds for q = 2 to a 
Greenberger-Horne-Zeilinger (GHZ) state 

I^gHI+^ + I-)^, (8) 

where | — ) = |0) — |1), and also this result leads to the well 
known closed formula for the partition function. Finally, the 
2D rectangular lattice with open (periodic) boundary condi- 
tions leads for q — 2 to a planar (toric) code state as intro- 
duced in the context of topological quantum error-correction 
fTjll . We remark that the connection between the 2D Ising 
model and planar (toric) code states was first proven and uti- 
lized in Ref. Q. 

3. — Solutions using quantum techniques. Having es- 
tablished the above formulation of the partition function Z G 
in a quantum language, we now show how to use powerful 
techniques established in quantum information theory to study 
Z G . Our general approach will be to translate insights regard- 
ing the states \4>g) into insights regarding the partition func- 
tions Z G . We remark that it is not obvious that the formula- 
tion in a quantum language simplifies the problem, in particu- 
lar as it involves quantum states of N q-dimensional systems, 
which are described by q N -dimensional vectors and thus have 
an increased number of parameters. However, we can make 
use of the following two advantageous features of the corre- 
spondence 0. First, as pointed out above, the states \4>g) 
belong to the class of stabilizer states, which are in general 
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a highly manageable family of multi-particle quantum states 
which has been extensively studied, and of which many prop- 
erties are known. Second, it follows from © that the state 
\4>g) is on ly used t° obtain an encoding of the interaction pat- 
tern of the model — i.e., the (oriented) graph G — and that all 
information regarding the specifics of the interactions — i.e., 
the Hamiltonians h e — is encoded into the single particle states 
\a e ). Translating properties of the states \4>g) into properties 
of the corresponding partition functions thus yields a way of 
unifying properties of models which are specializations of the 
current model (such as the Ising, Potts and clock models — see 
above), as one simply needs to consider different local states 
\a e ) for the same state \4>g)- 

Next we give three illustrations how to use stabilizer tech- 
niques to investigate the partition function Zq. 

3.1. Efficient evaluation of Zq. First we consider the prob- 
lem of evaluating Zq in a given set of hamiltonians {h e }. This 
problem is known to be NP-hard on general graphs. Efficient 
algorithms for special instances of graphs do exist, e.g., for 
the Ising model on planar graphs |2|. Here we show that: 
The partition function Zq can be evaluated efficiently (i.e., 
in polynomial time in the number of edges N) on all graphs 
G of logarithmically bounded tree-width Moreover, we 
provide an explicit algorithm based on the efficient classical 
description of quantum states in terms of tree tensor networks. 

To obtain this result, we use that any quantum state can be 
represented by a tree tensor network (TTN) fl6J], i.e., as a col- 
lection of tensors of rank Xi that are arranged according to a 
tree structure. In general, the maximal rank \ : = max^ of 
the tensors in the network grows exponentially with the num- 
ber of systems N. As shown in Ref. |6J], using the TTN de- 
scription one can extract information from the quantum state 
with overhead 0(N -poly(x))- In particular, one can calculate 
overlaps with product states. For stabilizer states on qubits 
(q = 2), it was shown how to obtain the optimal TTN (i.e. 
with minimal y) i n Ref. |2B, and this method can be extended 
to arbitrary q I19B . It was shown in Ref. |7] that, for a given 
stabilizer state the optimal \ is obtained by computing 
the (exponential of the) Schmidt-rank width of More- 
over, the optimal TTN description of the state can explic- 
itly be computed in poly(N, x) steps on a classical computer. 
Evaluation of overlaps with product states, and, hence, spe- 
cializing to the case of |V>g)> evaluation of partition functions 
Zq corresponds to a simple contraction of tensors, which can 
be done in 0(N ■ x 3 ) steps |0]. For \ipG)> one shows that the 
optimal x scales as the exponential of the tree-width t(G) of 
the graph G, by using that the Schmidt-rank width of \4>g) is 
proportional to the branch-width b(G) fl5ll of the underlying 
code Cq(o-, q) | 19]; further, one uses the inequalities 

fi(b(G)) < t{G) < h(b(G)), (9) 

where fi, are linear functions 1 15]. It follows that when- 
ever t(G) grows at most logarithmically with N, then x scales 
polynomially with N, and hence the partition function Zq can 
be evaluated efficiently. 



Note that the computationally hardest part in the above 
algorithm is to obtain the TTN description of \i/jq), while 
the calculation of overlaps with arbitrary product states, and 
hence the variation of coupling strengthes or different models 
for a fixed geometry, only scales linearly with N. We also 
emphasize that nonplanar graphs (of logarithmically bounded 
tree-width), as well as non-local interactions, are covered by 
this result. Results regarding efficient computation of homo- 
geneous Potts model partition functions on graphs of (loga- 
rithmically) bounded tree-width have been obtained before 
(l(l, though with entirely different methods. We emphasize 
that our approach, in contrast to previous approaches, can han- 
dle without difficulty also inhomogeneous models. Moreover, 
our method leaves a lot of space for generalizations. 

3.2. Dualities for planar graphs. Next we show how the 
partition function Zq of a planar graph G can be related to the 
partition function Zjj of its dual graph D Hcl. First, note that 
every orientation a of G induces an orientation of its dual D, 
which we also denote by a. We refer to Ref. flflll (p. 168) for 
details. Letting B(G a ) and B(D a ) be the incidence matrices 
of G a and D a , respectively, one then has B(G a )B(D' J ) T = 
1 10] (p. 169), which implies that the spaces Cq(ct, q) and 
C>D{p,q) are each other's duals, i.e., Cq(<7, q) 1 - — Cd{ct, q). 
Now, let F be the quantum Fourier transform, or generalized 
Hadamard operation, 

F:=^- J2 ^ kj \j)(k\, (10) 
V 9 j,k=o 

defined on C 9 . This operator satisfies FXF^ = Z and 
FZF^ = X, where Z and X are defined in ©. Using this 
property and the duality of Cg(ct, q) and Cn(cr, q), one finds 
that 




In other words, the states \iPg) an d \iPd) ore equal up to 
a Hadamard operation applied simultaneously on all local 
Hilbert spaces in the system. This result immediately relates 
the partition functions Zq and Zq, as one finds that 




where now \a' e ) := W\a e ). This leads to Zc(q, a, {h e }) — 
Zp(q, a, {h' e }), where the h' e are defined through 

e -PKU) :=-Ly e-l^e^W (13) 

for every j = 0, . . . , q — 1. 

As an example, consider the Potts model on the graph G, 
which is obtained by putting h e (0) := —J e and h e (j) = 
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whenever j ^ 0, for every e G E. The identity Jl 3i then 
yields 



g/3Je + q _ I if j = 



if j = l,...,q- 1. 



Defining interaction strengths J' e through 

e 0Je +q -l 



=/3J e 



1 



(14) 



(15) 




or, equivalently, (eP J <> — l)(e^ Je — 1) = q, one recovers the 
well know high-low temperature duality relation for the Potts 
model partition function 1 1], where the partition function on 
the graph G with interaction strengths { J e } is related to the 
partition function on the dual graph D with interactions { J' e }. 

3.3. Symmetries. Finally, general symmetries can easily 
be obtained. Here the essential observation is the following: 
any local unitary operator U := &) eeE U e having \4>q) as an 
eigenstate with eigenvalue A, yields a symmetry 



(16) 



where \a' e ) := XU e \a e ). A symmetry is obtained in the 
sense that different configurations of interactions (which are 
encoded in the Joy)) on the same graph G yield the same par- 
tition functionl 17]. Here again the fact that the states \4>g) are 
stabilizer states is particularly advantageous, as every \ipo) is 
the joint eigenstate of the q N local Pauli operators X(u)Z(v) 
(see © and below), which constitute the stabilizer group of 
IV'g)- Thus, this stabilizer group corresponds to a group of 
q N symmetries of Zq. 

Let us give an example of such a symmetry for q — 2. Let 
G be an arbitrary graph, let a be one of its vertices, and let 
E(a) be the set of edges incident with a. Let X[a] be the 
TV-qubit correlation operator on such that the the eth tensor 
factor is equal to the Pauli matrix X if e G E(a), and equal 
to / otherwise. Then X[a] stabilizes \tpG,2) and, using the 
relation with Zq, one obtains a relation between Zq(2; { J e }) 
and Z G {2; {J' e }), where J' e := - J e if e G E{a), and J' e := J e 
otherwise. 

4. — Relation to measurement based quantum compu- 
tation. So far we have used stabilizer techniques to study 
the partition function Zq, using the relation 0. Here we 
remark that this connection can also naturally be invoked to 
gain insights in the opposite direction — i.e., insights in Zq 
can be used to understand the properties of the states \^>q) 
and their possible role in quantum information tasks. In par- 
ticular, one may consider the question whether the states \iPg) 
are useful resources for measurement based quantum compu- 
tation (MQC), similar to the 2D cluster state, which in fact al- 
lows universal computation by local measurements only 1 18]. 
Like all stabilizer states, the \ipG) are natural candidates for 



resources for MQC, as these states are generally highly entan- 
gled and as they can be produced efficiently by a poly-sized 
quantum circuit. The connection with the partition function 
immediately relates the difficulty of evaluating Zq with the 
difficulty of classically simulating local measurements on the 
state \iPq). In particular, one finds that all (classes of) graphs 
for which the computation of Zq is hard, give rise to resource 
states for which simulation of MQC is a hard problem. For ex- 
ample, nonplanar graphs of large tree-width would typically 
give rise to states where MQC is NP-hard to simulate clas- 
sically. Conversely, known results regarding efficient evalua- 
tion of Zq on certain graphs — e.g., planar graphs for the Ising 
model (q — 2) — may lead to states \ipG) on which MQC can 
be simulated efficiently classically. We refer to Ref. fl4Tl . 
where this approach was first adopted to show that MQC on 
Kitaev's planar code states can be simulated efficiently classi- 
cally. 

5. — Summary and outlook. We have related the evalu- 
ation of the partition function on a general class of classical 
spin models to the calculation of overlaps between quantum 
stabilizer states encoding the interaction pattern, and product 
states encoding the details of the interaction. This allowed us 
to use powerful techniques established in quantum informa- 
tion theory to obtain novel classical algorithms to evaluate the 
partition function, and to investigate general symmetries and 
duality relations. We emphasize that in all cases we consider, 
the reformulation in terms of a quantum problem is virtual, in 
the sense that we provide an alternative way to obtain a clas- 
sical solution of the initial classical problem by making use 
of techniques established to describe and simulate quantum 
systems by classical means. 

Finally, we note that other classical combinatorial problems 
can be formulated in the language of quantum physics using a 
similar construction as outlined above. As an example, we 
mention that the weight enumerator of a classical code — a 
central quantity in classical coding theory — can be evaluated 
with similar methods |19]. Further, as the Potts model par- 
tition function is intimately related to problems in graph and 
knot theory through the relation with the Tutte polynomial, 
the present results also have implications in these fields. 
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